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Abstract 

We investigate contributions of Kaluza-Klein (KK) graviton in extra dimension models 
to the process pp —)■ pyp —)■ pyjX, where a proton emits a quasireal photon and is 
detected by using the very forward detectors planned at the LHC. In addition to the 
'yq initial state as in the Compton scattering in the standard model, the jg scattering 
contributes through the t-channel exchange of KK gravitons. Taking account of pileup 
contributions to the background and examining viable kinematical cuts, constraints on 
the parameter space of both the ADD (Arkani-Hamed, Dimopoulos and Dvali) model 
and the RS (Randall and Sundrum) model are studied. With 200 fb“^ data at a center- 
of-mass energy of 14 TeV, the expected lower bound on the cutoff scale for the ADD 
model is 6.3 TeV at 95% conhdence level, while a lower limit of 2.0 (0.5) TeV is set on 
the mass of the hrst excited graviton with the coupling parameter k/Mpi = 0.1 (0.01) 
for the RS model. 



1 Introduction 


After the discovery of the Higgs boson [y|^, the most important task of high energy physics 
at the LHC is to look for signatures of new physics beyond the standard model (SM). A 
possibility of introducing spatial extra dimensions has been discussed to explain a large 
hierarchy between the Planck scale (~ 10^® GeV) and the Fermi scale (~ 10^ GeV). Two 
representative scenarios in extra dimensional models are the large extra dimension model by 
Arkani-Hamed, Dimopoulos and Dvali (ADD) and the warped extra dimension model 
by Randall and Sundrum (RS) [^. The main difference between the two models is the 
role of extra dimensions solving the hierarchy problem while a common phenomenological 
consequence of the two models is Kaluza-Klein (KK) excitations of graviton. Therefore 
searching for KK gravitons in high-energy collider experiments is crucial to probe if our 
spacetime is four-dimension or more. 

At the LHG, signatures of ADD gravitons have been often studied in monojet or 
monophoton plus missing energy hnal states, where the missing energy is carried by 
real emissions of the KK gravitons. Moreover, ADD signatures from virtual KK graviton 
exchange have been sought in diphoton 10,11 , dilepton 12,13 and dijet 14 hnal states. 
For the RS model, on the other hand, the resonant production of the hrst excitation of 


the KK gravitons has been looked for in diphoton 10,11 , dilepton 15,16 , dijet 17 and 


di-VF boson 18,19 hnal states. No signihcant excess over the SM background has been 
observed so far, setting the limit, e.g. on the Planck scale in 3 -|- 5 spatial dimensions (M^) 
at around 4 TeV with 5 = 4 Im and on the scale of virtual graviton exchange at around 


7 TeV 14 for the ADD model, and on the hrst KK graviton mass at around 2.7 TeV with 


the coupling parameter k/M-pi = 0.1 for the RS model [1^. Thus an expectation to hnd the 
signal of KK graviton is postponed until the next (or future) stage of LHG experiments such 
as pp collision in y/s = 13 TeV or at the high-luminosity LHG with an integrated luminosity 
Lint = 3000 fb“^. Among the upgrade plans of LHG experiments, the installation of very 
forward detectors in addition to the central detectors to study the pp scattering with very 
large pseudorapidity of protons is proposed by ATLAS and GMS 20 . The former is called 


ATLAS Forward Physics (AFP) project, and the latter is called GMS-TOTEM forward 
detector scenario (GMS-TOTEM, in short). The forward detectors will be set at 220 m 
and 420 m from the interaction point in the AFP project while it is set at 420 m from the 
interaction point in addition to the TOTEM detectors in the GMS-TOTEM scenario. These 
forward detectors will detect the scattered protons which have momentum fraction loss 
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where p and p' denote the three momenta of incoming and outgoing protons, respectively. 
The acceptance of ^ is aimed as 0.0015 < .^ < 0.15 for AFP and 0.0015 < ^ < 0.5 for 
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Figure 1: Schematic diagram for pp —>■ p'yp —)• p^jX. 

CMS-TOTEM. 

In this paper, we study the possibility of using the forward detectors to hud indirect 
constraints on the KK graviton in both the ADD and RS models through the t-channel KK 
graviton exchange effects in 


pp —)■ P'yp —)■ P'jjX. 


( 2 ) 


The schematic diagram is shown in Fig. In this process, a quasireal photon 7 with low 
virtuality is emitted from a proton which scatters with a parton in the proton coming from 
the opposite direction. The proton which emits 7 is scattered with very small angle without 
dissociation into partons, and detected by the forward detectors. An earlier work to study 
the KK graviton in the above process has been presented in Ref. 21 , where the authors 


obtained the bounds on KK graviton from yg — yg (g is a quark or antiquark in a proton) 
processes performing a simple parton-level analysis. It must be noted that the process 
yg —)■ yg is the only subprocess of the process ([^ at the leading order in the SM. However, 
the exchange of the KK gravitons allows yg —)• yg as another subprocess. In our study, 
therefore, we revisit the KK graviton contributions to the process (|^ including all possible 
subprocesses. Moreover, we take into account the SM background coming from pileup events, 
whose contribution can be signihcant in the future higher-energy run, pointed out recently in 
We study constraints on the parameter space of ADD and RS models by examining 


Ref. 22 


viable kinematical cuts and taking account of parton-shower and hadronization effects in the 
hnal states as well as detector effects to hnd more realistic bounds on the models, which 


have not been discussed in Ref. 21 


This paper is organized as follows. In Sec. we give a brief review of the ADD and RS 
models to hx our notation. The results of numerical analysis are given in Sec. Section 
is devoted to summary of our Endings. 
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2 Models 


In this section we briefly review the ADD and RS models, focusing on interactions of 
KK gravitons to the SM fields. 

In the both extra dimension models, there appear KK towers of gravitons, and the 


effective interaction Lagrangian with the SM fields is given by 23,24 


n 


( 3 ) 


where is the energy-momentum tensor of the SM fields (see, e.g.. Ref. 25 for the explicit 
forms), is the nth KK mode of the graviton, and A is the relevant interaction scale. 

The ADD model is given in the (4-|-5)-dimensional spacetime where the extra 5-dimensions 
are compactifled on a torus with a common radius R. The four-dimensional Planck scale 
Mpi is related to the (4 -|- 5)-dimensional fundamental scale Md as Mpj owing 

to the Gauss’s law. All the SM fields are expected to be confined on the four-dimensional 
spacetime and only the gravitational interaction propagates into the extra dimensions. The 
scale of the interaction in Eq. ([^ is given by 


A = Mpi = Mpi/y^ 2.4 X 10^® GeV, 


( 4 ) 


where Mpi is the reduced four-dimensional Planck scale. Setting the fundamental scale Md 
to be 1 TeV, two extra dimensions 5 = 2 imply a large size of radius R ~ 0(0.1 mm). After 
compactiflcation of extra dimensions, there are KK excitations of gravitons whose spacing 
is given by ~ 1/R. This leads to infinity in the KK graviton propagator after summing up 
all KK modes and we replace the propagator by the cut-off parameter A^ (of order Md) for 


5 > 2 by 23 
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where is the mass of the nth KK mode of the graviton]^ 

The RS model is a five-dimensional model where one warped spatial dimension y is 
compactifled on /Z 2 orbifold. The metric is given by 


ds^ = e — dy'^, 


( 6 ) 


where (h) ^ = 0,1,2,3) and k denote the Minkowski metric and the AdSs curvature, 
respectively. The model has two D 3 branes at y = 0 and ttPc, the former is called the 
Planck brane and the latter is called the TeV brane. The SM fields are confined in the TeV 


^This cutoff scheme is the so-called Giudice-Rattazzi-Wells convention 23 , while there are other conven¬ 


tions, e.g. by Han-Lykken-Zhang 24 and Hewett 26 . See Ref. 13 for more details. 
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brane and only graviton propagates into the fifth dimension]^ With this setup, the hierarchy 
between the Planck (Mpi) and the Fermi {Mw) scales is explained when kvc — 12, and the 
scale of the interaction in Eq. ([^ is 

A = (7) 

Therefore the interactions of all the KK gravitons to the SM fields are suppressed by A-j^ ~ 


(P(TeV). The mass of the nth KK mode of the graviton is given by 27 


( 8 ) 


where Xn is a root of the Bessel functions of the first kind. The denominator of the graviton 
propagator is normal: 

1 




(9) 


where T^ denotes the graviton decay width 24 For the parameter scan, we use 

(3 = k/M-pi and me = nii, (10) 

which are commonly chosen. 


3 Numerical Analysis 

3.1 Signals 

Based on the RS graviton implementation 29 , we implemented the propagator in Eq. (|^ for 


the ADD gravitons by modifying the UFO file 30 32 , while we introduced higher KK gravi¬ 
ton modes for the RS model into FeynRules (3^. We use MadGraph5_aMC@NLO |34 


to generate parton-level events both for the signal and the SM background by employing 
the MSTW2008 PDF 35 through the LHAPDF interface (3^, with the factorization scale 


fixed at five times the Z-boson massj^ We use the flux of quasireal photons emitted from 


a proton in MadGraph5_aMC@NLO, implemented in a similar manner with PDFs 37 


and detail it in the appendix. The partonic events are passed to PythiaS 38 for parton 
shower and hadronization. 

We conduct analyses for the LHC at \/s = 14 TeV. As minimal event selections, we 
impose cuts on the transverse momentum pp and the pseudorapidity rj for the photon and 
the leading jet as 


pJ>40GeV, |r/^|<2.5, > 50 GeV, < 3.0, 


( 11 ) 


^There have been some variants of the RS model in which some of the SM fields are allowed to propagate 
into the bulk. But we do not discuss such possibilities furt her in this paper. 

^The widths are computed by the decay package 28 for each parameter point. The values are, e.g. 


rn 2 . 3.4 = 4.2,26,80,180 GeV for (/3, ma) = (0.05,1.2 TeV). 

^See Ref. 21 for a more detailed discussion on the scale choices. 
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Figure 2: Left: Photon transverse momentum distribution of the signal plus background for 
the ADD model with the cutoff scale At = 6 TeV in pp —)■ p'jp —)• p^jX at a/s = 14 TeV. 
The contributions of each subprocess are also shown. Right; The p/ distributions for At = 5, 
6 and 7 TeV. 


where jets are reconstructed by the anti-/cr algorithm 39 with a distance parameter of 0.6. 


The jet-energy resolution is set to 10%. To show our results, we assume the integrated 
luminosity Lint = 200 fb“^, and take into account the so-called survival probability S = 0.7 
that the proton remains intact and the detection efficiency of the photon of = 0.8. 
Therefore, the number of events is N = a x L^t x S x e^. 

In Fig. [^we show the pt distribution of the photon for the ADD model with the cutoff 
scale At = 6 TeV. Here, we generate events for the signal (blue lines) and the background 
(shaded) independently, and compare with the full sample (red solid) including the inter¬ 
ference between them. While only the yg initial state contributes in SM, the yg scattering 
contributes in extra dimension models, leading to about 25% enhancement of the ADD sig¬ 
nal rate. The signal dominates in the high-p'r region, and therefore a certain p/ cut largely 
reduces the background. We note that the interference between the signal and background 
is very small. In Fig. [fright) we show the p/ distributions of the signal plus background for 
At = 5, 6 and 7 TeV. 

Similar to the ADD case, the high-p^ photons are expected in the RS case, shown in 
Fig.|^left), where a benchmark (/, me) = (0.05,1.2 TeV) is taken. In this case, the n = 1—4 
KK graviton masses are 1.2, 2.2, 3.2 and 4.2 TeV, respectively, and the cutoff scale Aj^ is 
6.3 TeV. The contribution from the yp initial state enhances the signal by about 25%. In 
Fig. [fright) we also show the p/ distributions of the signal plus background for me = 0.9, 
1.2 and 1.5 TeV. 

In the following analyses, we impose the high-p-r photon selection cut 


p/ > 600 GeV 


( 12 ) 


to remove the background. 
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Figure 3: Same as Fig. but for the RS model with (5 = 0.05 and me = 1-2 TeV (left) and 
me = (0.9,1.2,1.5) TeV (right). 


3.2 Background from pileup events 


At the naive parton level, the photon-induced Compton 7 g —)■ 7 g process only contributes 
to the background. However, the qg —)■ 7 g and qq —)• 'jg processes, denoted as pp —)■ jjX, 
can contribute rather signihcantly to the background as seen below. The LHC experiments 
operate under a very high-luminosity condition such that the multiple proton-proton inter¬ 
actions (pileup) take place in the same bunch crossing. The average number of interactions 


per bunch crossing (/r) around 50 is expected in the Run-H operation 40 . The majority of 


the pileup events consists of elastic scattering, single and double diffractive scattering and 
inelastic scattering. Although the forward proton is absent in pp —)■ yjX, the pileup events 
may produce forward protons in the hnal state which overlap with the hard scattering events. 
Therefore, the following two SM processes are considered as the background in this study: 

1 . photon induced process 'yp —)■ ■yjX (referred to 7 p —>■ 7 - 1 - j), 

2. overlap between pp —)■ 'jjX and pileup events (referred to pp —)■ 7 -|- j -|- PU). 

In order to evaluate the effect of pileup events, minimum bias events were generated with 


PythiaS 38 . The average number of pileup events is assumed to be 50, and the overlap 


of pp —)■ jjX and pileup events are simulated by randomly taking multiple pileup events 
from the minimum bias sample for each pp —)• 'yjX event. The cross section of the second 
process is about 16 nb at y/s = 14 TeV, four orders of magnitude larger than the hrst one 
(about 3 pb). This is expected for the hard interactions involved in these processes as the 
gluon-initiated process is dominant compared to photon-initiated process at the LHC. With 
50 pileup events on average, there are always multiple protons in the forward region either 
from diffractive production or inside the proton remnant in case of inelastic scattering. 

One way to reduce the 7 + j + PU contribution is to use the fact that the forward proton 
and the particles in the main detector (photon and jet) are produced by different proton- 
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Figure 4: Distributions of xi^ 2 /^ for the overlapping events of pp —)■ jjX and pileup con¬ 
tribution (solid line), and the ■yp —)■ 7 j sample (hlled histogram), where p]. > 600 GeV and 
0.0015 < ^ < 0.5 are imposed. 


proton interactions. In order to investigate this correlation, the momentum fraction of the 
proton taken by the partons (xi and X 2 ) is calculated from the four momenta of the photon 
and the jet reconstructed in the main detector as 




+ 7 '^ ’^)- 


(13) 


Here, Ep is the energy of the proton beam. The variables Xi and X 2 correspond to the 
momentum fractions for the partons moving toward positive and negative z direction, re¬ 
spectively. 

For the j + q —)■ J + q process, since the photon, jet and the forward proton originate from 
the same interaction, one expects either xi/^ ~ 1 or X 2 /^ — 1, where ^ is dehned in Eq. ([^. 
On the other hand, if they come from different interactions, xi and X 2 have no correlation 
with Figure l^shows the distribution of the ratio Xi^ 2 /i for the two background processes 
after requiring the minimal selection cuts ( [II| ) and the > 600 GeV cut ( [I^ . The forward 
protons are required to have 0.0015 < ^ < 0.5. While there is a good correlation for the 
IP ^ 1 + j process which creates a narrow peak at one, the pp ^ 'y + j + PU sample shows 
a broader spectrum. Although there should not be a correlation between xi ^2 and ^ in the 
pp ^ ■y + j + PU sample, the distribution shows a peak around one. This is because the 
distribution is produced by taking the proton which gives the xi/^ or X 2 /^ value closest 
to unity if multiple protons are present in the forward region, and indeed there are several 
protons within this ^ range with /i = 50. We require the event to have a forward proton 
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Figure 5: The (left) and ^ (right) distributions for the 7 + j + PU and 7 + j samples with 
the 0.0015 < ^ < 0.5 cut plus the xi^ 2 li ratio cut, where > 600 GeV is imposed. For the 
^ distribution, only the 7 + j + PU sample is shown. 


^ range 

^ range cut 

Xi, 2 li ratio cut 

overall 

(0.0015, 0.5) 

1.0 

3.9 

3.9 

(0.0015, 0.15) 

1.2 

13 

16 

(0.1, 0.5) 

1.0 

4.5 

4.5 

(0.1, 0.15) 

4.2 

12 

48 


Table 1: Rejection factors for the 7 + j + PU by applying the ^ range cut and 0.9 < < 

1.1. The factors are calculated with respect to the 7 + j + PU sample with > 600 GeV. 


with either 


0.9 < xi/^ < 1.1 or 0.9 < X 2 /^ < 1.1. 


(14) 


Figurej^illustrates the effect of the 0 : 1 ^ 2 /'^ ratio cut in Eq. (14) inthep(}, (left) and (right) 
distributions, where the selection of the photon p'^ > 600 GeV is applied. Figure [fright) 
represents the ^ distribution of the protons in pileup events, which are produced by diffractive 
scattering (small ^ value) or by the hadronization into a proton from particles in the proton 
remnant in case of inelastic scattering (large ^ value). In Fig. we take the ^ range cut 
as 0.0015 < .^ < 0.5 as a representative case, while we show the rejection factors for the 
7 + j' + PU sample with the different .^-range cuts in Table With this Xi^ 2 /i selection, the 
7 + j + PU background is reduced by a factor of 4 — 50, depending on the ^ range cut, while 
the photon-induced 7 -|- j background mostly remains. 
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Figure 6: The ^ distributions of the signal plus background for the ADD (left) and RS (right) 
models in pp —)■ p'yjX at a/s = 14 TeV, including the pileup contribution to the background, 
where p/ > 600 GeV and 0.9 < < 1-1 are imposed. 


3.3 Kinematical distributions 

Including the pileup contribution as the SM background, in Fig. [^we show the ^ distributions 
for the signal plus background after imposing the high-p^ photon selection cut (12) and the 


Xi^ 2 /i ratio cut (14). For some lower in the ADD model and smaller me in the RS 
model, the signal rates now become visible over the SM background. Since the graviton 
signals distribute mostly in the ^ > 0.1 region, in the following analyses we impose the lower 
^ cut for the further background reduction as 


0.1 < ^ < 0.5 for the CMS-TOTEM, (15) 

0.1 < ^ < 0.15 for the AFP. (16) 


The pt distributions of the photon after all the event selections for the ADD and RS 
models are shown in Figs. [^and|^ respectively, where we take the above different ^ ranges. 
These illustrate that the SM background contribution is suppressed to a reasonable level. We 
note that, although we fixed the p/ selection cut at 600 GeV in this work, some optimizations 
may help to reduce the background further. 


3.4 Limits on the model parameters 

Finally, we would like to constrain the parameter space for each model. The expected exclu¬ 
sion limit on the parameter space in each model is derived by assuming a null observation. 
Taking into account only the statistical uncertainty, the function is dehned from the num¬ 
ber of signal and SM background events as = (Vs+s — NbY/Nb- The number of events 
are normalized to integrated luminosities from 20 fb“^ to 200 fb“^ for the ADD model, which 
has only one model parameter, i.e. the cutoff scale A^. For the RS model, on the other hand. 
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Figure 7: The distributions of the signal plus background for the ADD model with the 
different cutoff scales, = 5, 6 and 7 TeV, in pp —)■ pjjX at y/s = 14 TeV for the 
CMS-TOTEM (left) and the AFP (right), where p]. > 600 GeV and 0.9 < xi^ 2 /i < 1-1 are 
imposed. 




Figure 8: Same as Fig. but for the RS model with the hxed coupling parameter (3 = 0.05 
and the different KK graviton masses, me = 0.9, 1.2 and 1.5 TeV. 


we £x an integrated luminosity at 200 fb ^ and scan the two model parameters, i.e. f3 and 


me in Eq. (10). The model parameters are considered as excluded at 95% conhdence level 
(CL) if > 3.84 is satished. 

The is calculated for each point in the parameter space where the signal event was 
generated. The exclusion limit is determined at the parameter value crossing = 3.84, as¬ 
suming that the varies smoothly with the model parameters. The uncertainty of expected 
exclusion limit is evaluated by varying Ns+b and Nb within their Poisson uncertainties, and 
repeating the procedure to obtain the limit. With the kinematic selection described above, 
Nb = 13246 (1245) for 0.1 < .^ < 0.5 (0.1 < .^ < 0.15) is obtained for an integrated 
luminosity of 200 fb“^ which allows a reasonable estimate of the statistical uncertainty. 

The expected lower bound of A^^ for the ADD model is shown in Fig. as a function 
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Figure 9: The 95 % CL lower bound on the cutoff parameter in the ADD model for 
0.1 < .^ < 0.5 (red, CMS-TOTEM) and 0.1 < ^ < 0.15 (blue, AFP) as a function of the 
integrated luminosity. The shaded bands around the limit indicate the ±la uncertainties. 


of an integrated luminosity. With 200 fb“^ data, a lower bound of A = 6.3 (5.4) TeV can 
be achieved for the CMS-TOTEM (AFP). Selecting events in 0.1 < .^ < 0.5 gives a better 
results compared to a narrower range 0.1 < .^ < 0.15. These limits are comparable with the 


current dijet analysis 14 , where the s-channel virtual gravitons are considered. 


The exclusion limit for the RS model is shown in Fig. where an integrated luminosity 
of 200 fb“^ is assumed. Again, a more stringent limit is obtained by using the ^ range cut of 
0.1 < ^ < 0.5 (CMS-TOTEM). A 95% CL lower limit on the hrst KK graviton mass for the 
CMS-TOTEM is expected to be 2.0 (0.5) TeV for the coupling parameter = 0.1 (0.01), 
which is slightly weaker than the current limit from the dilepton analysis 16 . The limit for 
the AFP is weaker and the values are 1.7 (0.4) TeV for [3 = 0.1 (0.01). 


4 Summary 


We have studied possibilities to search for the KK graviton in the process pp —>■ pyp —)■ 
prfjX using the very forward detectors planned at the LHC. This process consists of two 
subprocesses yg —)• yg and yg —)• yg, where the latter appears only through the t-channel 
KK graviton exchange and has been overlooked in the previous study 21 . We examined 


all possible subprocesses and viable kinematical cuts to reduce the SM background to the 
signal events. 

A serious background to the signal process besides the SM Compton process comes from 
the overlap between hard scatterings, gg —)■ yg and gg —?• yg, and forward protons from the 
pileup events. We showed that such background events could be reduced by requiring a ratio 
of Xi/^ (or X 2 /i) to be between 0.9 and 1.1. 

Taking account of parton-shower and hadronization effects in the hnal state, we found 
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Figure 10; The 95 % CL lower bound on the hrst KK graviton mass me as a function of the 
coupling parameter = k/Mpi in the RS model for an integrated luminosity of 200 fb“^. 
The results are shown for 0.1 < ^ < 0.5 (red, CMS-TOTEM) and 0.1 < ^ < 0.15 (blue, 
AFP). The shaded bands around the limit indicate the ±1 ct uncertainties. 

realistic constraints on the model parameter space of ADD and RS models for the 14-TeV 
LHC. The 95% CL lower bound on the cutoff scale A^^ in the ADD model is expected to 
be 6.3 TeV with the 200 fb“^ data for the CMS-TOTEM (5.4 TeV for the AFP). For the 
RS model, the lower bound of the hrst KK-graviton mass me is 2.0 TeV (0.5 TeV) with 
k/Mpi = 0.1 (0.01) for the CMS-TOTEM. Those for the AFP is 1.7 TeV (0.4 TeV) with 

k/Mpi = 0.1 ( 0 . 01 ). 

The process we studied in this article is not the conventional s-channel KK graviton 
productions, and hence the very forward detectors at the LHC might give us a new and 
complementary opportunity to search for new physics beyond the SM. 

Acknowledgements 

K. M. would like to thank the KEK theory group and Ochanomizu University for the warm 
hospitality, where a hnal part of this work was done. K. M. is supported in part by the 
Belgian Federal Science Policy Office through the Interuniversity Attraction Pole P7/37, 
and by the Strategic Research Program “High Energy Physics” and the Research Council of 
the Vrije Universiteit Brussel. 


Appendix: Photon flux from a proton 

The quasireal photon flux from a proton is described by the equivalent photon approximation 


(EPA) as 41 


/(e,Q 


' ) 
max/ 


OEM (1 — 0 


vr 




^ max 


Qm\n 


(17) 
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Figure 11: The EPA flux from a proton as a function of where the analytic expression 
in Eq. ( [ItI ) with Q^ax = 2 GeV^ (black) is compared with the numerical outcomes from 
MadGraph 5^MC@NLO for = 2 GeV^ (red) and M| GeV^ (blue). The ratios 
of the MadGraph5_aMC@NLO prediction to the analytic result are also shown in the 
bottom panel. 


where ^ = E^/Ep and Q^ax energy fraction of the photon and the maximal value of 

the integration over the virtuality of the photon g^, respectively. Here, g^in = HX 
and Ql = 0.71 GeV^. We usually choose g^ax ~ 2 GeV^ since the contribution above a few 
GeV^ is negligible. The function (p{x,y) is defined by 


^{x,y) = (1 + az) 


+ 


c(4 + z) 


In 


In 


1 + y 


y 

1 + y — b 


+ 5 : 


^ k{l + y)^ 
3 ik 


+ 


(1 — b)z 
Ay{l + yY 




1 + 1 / 2 ^^k{l + y) 


(18) 


with z = /(1 — x). The parameters a, 6, c are related to the mass and the magnetic moment 

of the proton and the values are given as a = 7.16, b = —3.96 and c = 0.028. 

"" "" ■ in a 


The above structure function is implemented in MadGraph5_aMC@NLO 37 


similar manner with PDF. In Fig. [T^ we show the photon flux as a function of comparing 
the above analytic function with the numerical outcome of MadGraph5_aMC@NLO. 
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